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I.  INTRODUCTION 


The  purpose  of  this  report  Is  to  present  a  mathemetlcal  model  of  the 
geometry  of  a  cylindrical  port  motor  cast  with  a  misaligned  mar.drel.  This  model 
was  developed  to  determine  the  burning  surface  area  and  free  volume  of  such 
motors. 

This  report  also  Includes  a  detailed  description  of  the  geometry  model.  In 
formulating  this  model,  two  basic  types  of  mandrel  misalignment  were  considered: 
mandrel  displacement  and  mandrel  cocking.  In  addition  to  the  model  description, 
two  appendices  are  included.  Appendix  A  presents  an  HP-4 1C  calculator  program 
and  Appendix  B  presents  an  example  of  the  application  of  the  geometric  model. 

The  details  presented  In  this  report  are  the  result  of  work  conducted  at  the 
Propulsion  Directorate  of  the  US  Army  Missile  Command.  The  purpose  of  this  work 
was  to  obtain  a  better  Insight  Into  the  geometrical  nature  of  mandrel 
misalignment. 

II.  GENERAL 

The  cylindrical  port  grain  Is  one  of  the  most  versatile  and  widely  used 
solid  rocket  motor  configurations.  This  motor  geometry  Is  widely  employed 
throughout  the  Industry.  One  of  the  more  common  applications  of  cylindrical 
port  grains  Is  In  subscale  ballistic  test  motors.  The  characterization  of  pro¬ 
pellant  burning  rates  Is  one  of  the  primary  uses  of  the  subscale  test  motors. 
Typically,  when  a  cylindrical  port  motor  Is  employed  In  burning  rate  charac¬ 
terization,  the  motor  Is  designed  with  a  burning  surface  area  profile  that  is 
essentially  constant  with  respect  to  web  distance  burned.  Thus,  when  fired,  the 
motor  will  operate  at  a  relatively  constant  chamber  pressure.  The  burning  rate 
of  the  propellant,  at  the  average  operating  pressure  of  the  motor.  Is  determined 
by  dividing  the  web  thickness  of  the  motor  by  the  burn  time.  This  entire  analy¬ 
sis  method  Is  based  on  the  assumption  that  the  web  thickness  of  the  motor  Is  a 
known  quantity.  Therefore,  It  Is  essential  to  this  method  that  the  web  distance 
be  uniform  over  the  entire  length  of  the  grain.  As  a  result  of  this  assumption, 
a  major  source  of  experimental  error  In  the  determination  of  burning  rate  from 
ballistic  test  motor  firings  Is  ballistic  test  motors  that  do  not  have  a  uniform 
web. 


The  major  cause  of  variations  In  the  web  thickness  for  cylindrical  port 
motors  Is  mandrel  misalignment.  Mandrel  misalignment  essentially  means  that 
when  the  motor  was  cast  the  axis  of  symmetry  of  the  mandrel  (and  thus  of  the 
motor  port)  did  not  coincide  with  the  axis  of  symmetry  of  the  motor  case.  This 
condition  causes  a  variation  of  the  web  thickness  over  the  length  of  the  grain 
which  means  that  the  burning  surface  will  not  contact  the  motor  case  wall  uni¬ 
formly.  As  a  result,  the  burning  rate  analysis  method  which  Is  based  on  the 
assumption  that  the  entire  burning  surface  contacts  the  motor  case  wall  at  the 
same  Instant  and  Is  rendered  useless. 
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Since  the  cylindrical  port  motor  is  such  a  basic  propellant  development 
tool,  it  is  essential  to  obtain  a  better  understanding  of  the  Influences  of 
mandrel  misalignment  on  the  performance  of  such  motors.  The  first  step  in 
obtaining  this  understanding  is  to  acquire  a  knowledge  of  the  geometry  of  misa¬ 
ligned  motors.  It  should  be  noted  that  the  effects  of  mandrel  misalignment  on 
the  performance  of  solid  rockets  were  extensively  investigated  by  Maykut  [1]. 

The  purpose  of  these  studies  was  to  investigate  the  effect  of  various  grain 
asymmetries  on  the  delivered  Impulse  of  a  rocket  motor.  In  these  studies  a 
generalized  grain  geometry  computer  code  was  employed.  One  feature  of  this  code 
was  the  ability  to  solve  for  the  surface  histories  of  various  asymmetric  pro¬ 
pellant  grains  [2].  While  this  code  was  capable  of  analyzing  the  geometry  of  a 
misaligned  cylindrical  port  motor,  the  general  nature  of  the  code  made  it 
somewhat  cumbersome  to  use.  As  a  result,  it  was  considered  advantageous  to 
Independently  develop  a  geometry  model  for  the  specific  class  of  motors  con¬ 
sidered  in  this  study. 

III.  MANDREL  MISALIGNMENT 

The  first  step  in  considering  mandrel  misalignment  in  a  cylindrical  port 
rocket  motor  is  to  consider  the  general  geometry  of  the  motor.  In  a  perfectly 
aligned  motor,  the  port  of  the  grain  and  the  motor  case  will  have  the  same  axis 
of  symmetry.  Figure  1  shows  the  geometry  of  such  a  motor.  The  problem  created 
by  mandrel  misalignment  is  that  the  motor  port  and  motor  case  do  not  have  a 
common  axis  of  symmetry.  In  order  to  begin  evaluation  of  the  nature  of  mandrel 
misalignment,  first  consider  the  case  where  the  port  and  case  axes  are  parallel 
but  do  not  coincide.  A  cross-section  of  the  motor  taken  through  a  plane  perpen¬ 
dicular  to  the  axes  of  symmetry  will  reveal  circular  port  and  motor  case  cross- 
sections.  These  circles  are  not,  however,  concentric.  As  the  propellant  port 
burns  out  radially  the  radius  of  the  port  will  Increase.  Eventlally  one  point 
on  the  port  cross-section  will  contact  the  case  wall.  This  point  defines  the 
region  where  the  misaligned  motor  differs  from  the  perfectly  aligned  motor. 

Until  the  point  of  contact  the  aligned  and  misaligned  motors  will  exhibit  the 
same  burning  surface  area  history. 

For  the  aligned  motor,  wall  contact  occurs  along  the  entire  periphery  and 
thus  indicates  the  time  of  motor  burnout,  while  for  the  misaligned  motor, 
wall  contact  creates  a  sliver  zone,  this  sliver  zone  is  the  cross-sectional 
area  of  propellant  remaining  at  the  point  of  first  wall  contact.  The  misaligned 
motor  will  continue  to  operate  as  the  sliver  zone  burns  out.  This  sliver  zone 
has  a  surface  area  that  decreases  as  web  distance  burned  Increases.  The  sliver 
will  result  in  an  extended  motor  tall-off  on  the  pressure-time  trace  for  the 
misaligned  motor.  Figure  2  presents  the  burning  profile  for  a  misaligned  cross- 
section. 

The  next  step  is  to  develop  a  mathematical  model  of  the  misaligned  cross- 
section.  Consider  the  misaligned  port  for  the  propellant  grain  at  a  given 
cross-section: 

The  ridius  of  the  propellant  grain  is  given  by: 

R  (t)  -  R(0)  +  T  (1) 


2 


PORT  AXIS  OF 


CASE  AXIS  OF 


Figure  2.  Cross-section  burn  profile  for  C-P  grain  cast  with 
misaligned  mandrel. 


Where: 

R  ( x)  is  the  radius  of  the  grain 
R  (0)  Is  the  Initial  grain  radius 
and  X  Is  the  web  distance  burned. 

The  intersection  of  the  propellant  port  and  the  motor  case  Is  given  by  the 
coordinates; 

(Xi,  +  Yi) 

If 

R  (x)  <  Rc  -  AX 


4 


There  Is  no  Intersection 


R  ( x)  Rc  -  AX 

Xi  =  R^  (t)  -  Rp^  -  AX^ 


^1= 

Where 

R(;  Is  the  Inside  radius  of  the  motor  case 

M  Is  the  magnitude  of  the  mandrel  offset 

Xj  Is  the  X-coordlnate  of  the  Intersection  and 

Yj  Is  the  Y-coordlnate  of  the  Intersection. 

The  perimeter  of  the  burning  surface  of  propellant  at  a  given  cross- 
sectional  plane  Is  given  by: 

P  (  x)  =  7t  0R  ( x) 

180 

Where: 


R  ( x)  ^  Rc  -  AX 
01  =  360® 


Xj;  <  -  AX 


9,  -  360-  -  2  T.n-1 


Xi  =  -AX 
01  »  180® 


I 


r 


01  2  TAN  ^ 

(8) 

Where: 

P  ( t)  Is  the  perimeter  of  the  propellant. 

The  cross-sectional  area  of  propellant  at  a  given 
given  by: 

cross-sectional  plane  is 

Acr  (t)  =  ^  (R^c  ^2  -  (t)  01^  +  2  Ai 

360 

(9) 

Where 

If 

R  ( x)  _<  Rj.  -  AX 

02  =  360* 

(10) 

Aj  =*  0 

(11) 

If 

Xi  <  0 

Y 

02  =*  360*  -  2  TAN"1 

(12) 

If 

X 

M 

H 

O 

02  *  180* 

(13) 

If 

Xi  >  0 

02-2  TAN"1  ll 

Xl 

(U) 

And 

Ai  -  [S  (S-AX)  (S-R(t))  (S-Rc)]^/2 

(15) 

Where; 

S  =  1/2  (AK  +  R  (t)  +  Rc)  (16 

Where: 

Acr  ( 1-s  the  propellant  cross-sectional  area 

At  a  cross-sectional  plane,  the  distance  for  the  shortest  propellant 
web  Is  given  by: 

T  sw  =  Rc  -  R(o)  -  ^  (^^^ 

The  web  distance  for  total  propellant  burnout  at  a  cross-section  Is 

given  by; 

t  pbo  =  Rc  -  R(0)  +  AX  (18) 

A  complete  cross-sectional  view  of  the  propellant  grain  Is  shown  In 

Figure  3. 


Figure  3.  Cross-sectional  view  of  C-P  grain  cast  with  an  offset  mandrel. 


IV.  MOTOR  GEOMETRY 


With  the  cross-sectional  geometry  of  the  propellant  grain  completely 
detailed,  the  next  step  is  to  consider  the  geometry  of  the  entire  motor.  In 
order  to  consider  the  motor  geometry  a  set  of  coordinate  systems  must  be 
established.  Two  coordinate  systems  will  be  considered,  one  for  the  motor  case 
and  one  for  the  mandrel.  Descriptions  of  the  coordinate  systems  are  as  follows: 

For  the  motor  case  - 

X  -  An  axis  in  a  plane  perpendicular  to  the  axis  of  symmetry  of  the 
motor  case 

Y  -  An  axis  in  the  same  plan  as  the  X-axis  and  perpendicular  to  the  X- 

axls  and  the  axis  of  symmetry  of  the  motor  case 

Z  -  The  axis  of  symmetry  of  the  motor  case. 

For  the  mandrel  - 

X  -  An  axis  in  a  plane  perpendicular  to  the  axis  of  symmetry  of  the 
mandrel 

Y  -  An  axis  in  the  same  plane  as  the  X-axis  and  perpendicular  to  the  X- 

axls  and  the  axis  of  symmetry  of  the  mandrel 

Z  -  The  axis  of  symmetry  of  the  mandrel. 

Two  possible  cases  of  mandrel  misalignment  will  be  considered.  The  first 
case  is  a  displaced  mandrel  and  the  second  is  a  cocked  mandrel.  The  following 
are  descriptions  of  the  two  resulting  motor  geometries. 

A.  Displaced  Mandrel 

In  the  case  of  the  displaced  mandrel,  the  assumption  is  made  that  the 
sides  of  the  mandrel  are  parallel  to  walls  of  the  motor  case  but  the  axis  of 
symmetry  of  the  mandrel  (Z-axls)  is  displaced  a  dlstanc^e  AX  from  the  axis  of 
symme^try  of  the  motor_case  (Z-axis).  Thus,  the  X  and  X  axes  are  colinear,  the  Y 
and  Y,  and  the  Z  and  Z  axes,  respectively,  are  parallel.  The  geometry  is  pre¬ 
sented  in  Figure  A. 
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For  the  displaced  mandrel  the  propellant  cross-section  at  each  Z- 
coordinate  is  the  same.  Therefore,  for  a  given  web  distance  burned  the  pro¬ 
pellant  perimeter  and  cross-sectional  area  are  constant  with  respect  to  Z. 

Thus,  the  propellant  burning  surface  area  is  given  by; 

Ab  (x)  =  L  (x)  P  (t)  +  Acr  (x)  Ngb  (19) 

Where 

R  (t)  is  given  by  Equation  (1) 

and  L  (x)  =  L(0)  -  2x  Neb  (20) 

Where 

Ab  ( x)  is  the  burning  surface  area  of  the  motor 
L  ( x)  is  the  length  of  the  grain 
L  (O)  is  the  initial  length  of  the  grain  and 
Neb  Is  the  number  of  ends  that  are  burning. 

The  free  volume  of  the  motor  is  given  by: 

V(x)  -  I'(O)  -  L  (x)  Acr  (x)  (21) 
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Where 


V( t)  is  the  free  volume  of  the  motor. 

B.  Cocked  Mandrel 

In  the  case  of  the  cocked  mandrel  two  general  geometries  will  be  con¬ 
sidered.  These  are  a  mandrel  that  is  cocked  at  the  top  of  the  motor  case  and  a 
mandrel  that  is  cocked  at  both  the  top  and  the  bottom  of  the  motor  case.  The 
following  presents  the  details  of  the  two  geometries. 

1.  Mandrel  Cocked  With  Respect  to  the  Motor  Case  Top 

In  the  case  of  the  cocked  mandrel  the  assumption  is  made  that  the 
axis  of  symmetry  of  the  mandrel  (Z-axls)  is  cocked  with  respect  to  the  axis  of 
symmetry  of  the  motor  case  (Z-axls).  In  the  case  where  the  mandrel  is  cocked 
with  the  respect  to  the  top  of  the  motor  case,  the  assumption  is  made  that  the 
coordinate  sy^t^m^  of  the  motor  case  and  the  mandrel  have  the  same  origin. 
However,  the  X-Y-Z  coordinate  system  is  created  by_rotating  the  X-Y-Z  system 
abou£  the  Y-axis.  Therefore,  the  X-,  Z-,  X-,  and  Z-axes  are  coplanar  and  the  Y- 
and  Y  -axes  are  identical.  The  geometry  is  presented  in  Figure  5. 


Figure  5.  Cocked  mandrel  configuration  (cocked  at  top). 
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In  order  to  determine  the  geometry  of  a  grain  created  with  a  cocked 
mandrel  three  simplifying  assumptions  are  implied.  These  are; 

a.  Axial  dlstances_along  the  propellant  grain  will  be  determined 
along  the  Z-axls  Instead  of  the  Z-axis. 

b.  The  propellant  cross-section  of  the  unburned  portion  in  the  X-Y 
plane  is  circular  Instead  of  elliptical. 

_  c.  The  propellant  burns  radially,  in  the  X-Y  plane  instead  of  the 

X-Y  plane. 


_  These  assumptions  are  justified  by  the  fact  tha^  the  angle  between 
the  Z  and  Z  axes  (which  is  the  same  angle  between  the  X  and  X  axes)  will  be  very 
small  and  thus  the  cosine  of  tlie  included  angle  will  be  very  close  to  unity.  In 
order  for  distances  along  the  Z-axls  to  exceed  distances  along  the  Z-axis  by 
more  than  .1%  the  Included  angle  must  exceed  2.5“ .  This  angle  should  be  well 
within  the  region  of  mandrel  misalignment  that  is  normally  encountered.  Thus, 
because  of  the  very  small  Included  angle  the  unburned  propellant  port  should  be 
essentially  circular  in  the  X-Y  plane.  Also,  this  small  included  angle  means 
that  web  distances  burned  along  the  X-axes  are  essentially  unchanged  when  pro¬ 
jected  on  the  X-axis.  And  finally,  the  effects  of  assumptions  a.  and  c.  above 
tend  to  cancel  each  other  and  thus  increase  the  accuracy. 


The  geometry  of  a  propellant  grain  cast  with  a  cocked  mandrel  can 
be  considered  to  experience  four  distinct  phases  as  the  motor  progresses  from 
the  initial  state  to  motor  burnout.  These  four  phases  are: 


PHASE  1.  The  port  of  the  propellant  is  totally  circular.  The 
short  propellant  web  had  not  burned  out  at  any  axial 
cross-section. 


PHASE  2.  The  short  propellant  web  has  burned  out  for  cross- 

sections  in  upper  portion  of  the  grain.  The  remainder  of 
the  grain  has  a  circular  port. 

PHASE  3.  The  short  propellant  web  had  burned  out  for  the  entire 
length  of  the  grain.  There  are  no  cross-sections  for 
which  total  propellant  burn  out  has  occurred. 


PHASE  4.  The  cross-section  at  the  bottom  of  the  motor  has 
experienced  total  propellant  burn  out. 


The  next  step  is  to  consider  the  geometry  of  the  motor  during  each  of  the 
following  four  phases: 
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PHASE  1 


0  £  X  £  Tl 

Where 

Rc  -  R(0)  -  AXx 


1  - 


Ntop^ 

L(0) 


AXt  -  AX  (Z-L(O)) 


(22) 


(23) 


Where 


Tj  Is  the  web  distance  burned  for  short  web  bum  out  at  the  top  of 
the  grain. 

Ntop  ®  0  If  the  top  end  Is  Inhibited 

*  1  If  the  top  end  Is  uninhibited 

and  AX-j-  Is  the  Initial  off-set  of  the  mandrel  axis  at  the  top  of  the  grain. 

The  burning  area  of  the  motor  Is: 

Ab(x)  =  P(t)  L(x)  +  (Nbot  +  Ntop>  Acr(t)  (24) 

Where 

L  (x)  =.  L  (0)  -  X  (Nbot  +  Ntop)  <25) 

Nbot  =  0  If  the  bottom  is  Inhibited  (26) 

=■  1  If  the  bottom  Is  uninhibited  (27) 

and  for  all  phases 

R( x)  is  determined  from  Equation  (1) 

The  free  volume  of  the  motor  Is  given  by  Equation  (21). 

PHASE  2 

xi  £  X  <  X2 

Where 

Rc  -  R(0) 


(28) 


1  + 


Nbot 
L(0) 
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^2  is  the  web  distance  burned  for  short  web  burn  out  at  the  bottom 
of  the  grain 

The  burning  surface  area  of  the  motor  is  given  by: 


Ab  ( x)  -  P  (x,  Zbot)  (^ub  “  ^bot)  ^  (T,z)dz 


+  Nbot  Acr  2bot^  "*■  ^top  Acr  ^top^ 


Where 


Zbot  “  i®  Z-coordlnate  of  the  bottom  of  the  grain 

Zub  “  is  the  Z-coordinate  at  which  the  cross-section  is  at  the 
exact  point  of  short  web  burn  out  and 


Ztop  ~  is  the  Z-coordinate  of  the  top  of  the  grain. 


Where 


Zbot  ( t)  =  X  Nbot 

Zub  (x)  =  L(0)  (Rc  -  R(0)  -  x) 


Ztop  ( ~  L(0)  X  Nfop 
and  note  that  for  all  phases: 


M(Z)  =  AXx 


The  Integral  term  can  be  approximated  by  applying  the  trapezoidal  rule 
over  11  points: 


“top 

f  ^ 

/  P(  x,z)  dz  =  —  /  ^ 


(p  (u.Zi)  +  P  (x,  Zi  -  AZ)) 


Where; 


AZ  = 


(T)  -  Z  XT) 

top _ UP 

10 


Zl  =  Zub  +  i  (AZ)  i  =  1,  2 . 10 

Thus,  the  burning  surface  area  is  given  by; 


(35] 

(36] 


10 

Ab<^)  ■  P(^.  Zbot)  (Z„b  -  Zbot)  .  “  E 

i=l 


+  Nbot  Acr  (t»  Zbot)  ^top  (t»  ^top) 

The  free  volume  of  the  motor  is  given  by; 

^top 

V(t)  =  L(0)  -  Acr  (t.  Zbot)  (Zub  “  ^^ot)  “  /  ^cr  (t,Z)  dz 

Zub 

This  can  be  approximated  by; 

V(  t)  =  uRc^  “  ^cr  (^»  ^bot)  (^ub  ^hot^ 

10 

_4|.  V  ('•  ^‘)  *  *er  t,  Zi  -  4Z)) 
i=l 

The  grain  length  is  given  in  Equation  (25). 

PHASE  3 


(37] 


(38) 


(39) 


T2  ^  T  <  T3 
Where 

Rc  -  R(0) 


T3 


1  - 


AXt  Ntop) 
L(0) 


Where 


(40) 


T3  is  the  web  distance  burned  for  total  propellant  burn  out  at  the 
bottom  cross-section. 


The  burning  surface  area  of  the  motor  is  given  by: 


^top 


This  can  be  approximated  by: 


lU 

AbCt)  °  T  ^  (pC-t.Zi)  +  P  (t.  Zi  -  Az)) 

i=l 

■*■  ^bot  Acr  Zbot^  ^top  Acr  (^pZtop) 


Where 


Zbot  ( determined  from  Equation  (30)  and 
Ztop  ( "t)  is  determined  from  Equation  (32) 

Ztop  (t)  -  Zbot  ^ 


and  Zi  =  Zbot  +  i  (AZ)  i  =  1,  2,  .  .  10 

The  free  volume  of  the  motor  is  given  by: 


V  (t)  =  tcRc2l(0)  -  /  A  (T,Z)  dz 


This  can  be  approximated  by: 


V(  t)  =  2  ^cr  ^T^jZi^  +  Aj,^.  ,  Zj^  -  Azjj 


The  grain  length  is  given  in  Equation  (25). 
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PHASE  4 


”^3  ^  ^  ^mbo 

where 


%bo 


Rco; 


/,  ,  Ntop  \ 

y  L  (0)  j 


(47) 


Where 

x^tjQ  is  the  web  distance  burned  for  total  motor  propellant  burn  out. 

The  relationships  for  burning  surface  area  and  motor  free  volume  are  the  same  as 
those  presented  for  Phase  3  with  the  exception  that; 

(48) 


(49) 

the  motor  case  bottom  and  top 

A  variation  of  the  cocked  mandrel  geometry  can  be  achieved  by  con¬ 
sidering  the  case  where  the  matidr^l  is  cocked  at  both  the  top  and  bottom  of  the 
motor  case.  In  this  case  the  Z-axis  is  created  by  rotating  the  Z-axis  about  an 
axis  which  is  parallel  to  the  Y-axis  and  that  passes  through  the  centroid  of  the 
unburned  propellant  grain.  This  geometry  is  shown  in  Figure  6. 

The  geometry  of  propellant  grain  can  be  determined  by  applying  the 
relationships  derived  from  the  situation  where  mandrel  is  cocked  about  the 
bottom  of  the  motor  case. 

The  burning  surface  area  of  the  propellant  grain  is  given  by; 

Ab  (  t)  =  2  Ajj  (t,  L(0)  ,  N{)ot^  (5Q) 


Zbot  (t)  =  (r  -  Rc  -  R(0)) 


L(0) 

AXt 


The  length  of  the  propellant  grain  is  given  by; 
L(  x)  =  Zfop  ( ''•)  ~  Zbot  ^ 

2.  Mandrel  cocked  with  respect  to  both 


Where; 

Ab  (x,  L(0),  Nbot) 

is  the  surface  area  determined  for  a  propellant  grain  created  by  a  mandrel 
cocked  at  the  top  only. 
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Figure  6.  Cocked  mandrel  configuration  (cocked  at  bottom  and  top) 


The  Inputs  to  the  burning  surface  area  relationships  are: 

l’(o) 

L(0)  -  -  (51) 

2 

and 

Nbot  =*  0  (52) 


Where 

L'(0)  Is  the  Initial  length  of  the  propellant  grain  created  by  a  mandrel 
cocked  at  both  the  bottom  and  top. 

Likewise  the  free  volume  of  motor  is  given  by: 

V(t)  -  2V  (t,  L(0),  Nbot)  (53) 

Note  that  these  relationships  apply  only  for  the  cases  where  either  the  top 
and  bottom  of  the  grain  are  both  inhibited  or  both  uninhibited. 


Thus 


If  Nfop  “  0  both  ends  are  inhibited  (54) 

if  Nj-qp  =  1  both  ends  are  uninhibited  (55) 

Also,  note  that  the  geometry  for  grains  which  were  generated  by  cocking  the 
mandrel  about  horizontal  axes  located  on  various  points  on  the  Z-axis  can  also 
be  determined  from  the  previous  relationships.  These  results  can  be  obtained 
by  adding  the  results  for  two  appropriate  motor  geometries  which  were  cocked  at 
the  top  of  the  motor  case. 

V.  CONCLUSIONS 

The  mathematical  model  presented  in  this  report  provides  a  means  to  deter¬ 
mine  the  geometrical  profile  of  cylindrical  port  motors  cast  with  misaligned 
mandrels.  This  model  should  serve  as  a  valuable  tool  in  determining  the  effect 
of  mandrel  misalignment  on  the  pressure-time  traces  of  ballistic  test  motors. 

In  this  application  the  model  could  be  used  to  make  some  determination  on  the 
accuracy  of  burning  rate  data  obtained  from  motors  with  various  degrees  of  misa¬ 
lignments.  Thus,  this  model  could  be  used  to  establish  a  set  of  criteria  for 
the  accuracy  of  burning  rate  data  obtained  from  cylindrical  port  ballistic  test 
motors. 
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APPENDIX  A 


HP-4 1C  PROGRAM 

The  mathematical  model  presented  in  this  report  has  been  incorporated  into 
a  program  for  an  HP-4 1C  calculator.  This  appendix  is  intended  to  provide  all 
the  information  required  to  Install  and  operate  this  program.  This  program, 
when  Installed  on  an  HP-41  calculator  system,  will  prove  to  be  a  useful  analysis 
tool.  The  program  as  presented  will  provide  the  user  with  a  convenient  and 
accurate  method  for  evaluating  the  geometry  of  misaligned  cylindrical  port 
motors.  The  following  provides  complete  ooerating  instructions,  a  set 
of  sample  problems,  and  a  listing  of  the  program.  Also  provided  is  all  the 
required  storage  register  and  calculator  status  information  needed  to  implement 
the  program. 

A.  Operating  Instruction 

In  order  to  Implement  the  program  presented  in  this  report  the  following 
equipment  is  required: 

1  -  HP-41CV  calculator 
or 

1  -  HP-4 1C  calculator  with  1  HP  82170A  quad  memory  module 

1  HP  83143A  thermal  printed/plotter 
or 

I  -  HP  83162A  thermal  printer/plotter  with  HP  82160A  HP-IL  module 

To  operate  the  program  the  printer  should  be  mated  with  the  calculator  in 
the  appropriate  manner.  The  calculator  should  then  be  configured  to  size  43  and 
placed  in  the  user  mode.  Table  A-1  provides  a  step  by  step  key  sequence 
required  to  operate  this  program. 
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TABLE  A-1.  Program  Instructions 


D 

I 

C 


i 


i 


STEP 

INSTRUCTIONS 

INPUT 

1. 

Load  Program. 

2. 

Clear  all  resisters. 

3. 

Initialize  program. 

4. 

5. 


5. a 


5.b 


S.b.l 


3.b.l .a 


5.b.l .b 


5.b.2 


5.b.2.a 


Key  in  case  radius,  R^. 

Indicate  if  1  or  0 

the  mandrel 
is  cocked. 

If  mandrel  is  1 

cocked. 

If  mandrel  is  0 

not  cocked 
go  to  Step  6. 

If  mandrel  is  1  or  0 


cocked  indicate 
if  it  is  cocked 
at  the  top  only. 

If  mandrel  is  1 

cocked  at  top  only. 

If  mandrel  is  0 

not  cocked  at 
the  top  only, 
go  to  Step  5.b.2, 

If  mandrel  is  N^ot 

cocked  at  top  only 
indicate  if  the 
bottom  is  burning. 

If  mandrel  is  Ntop 

cocked  Indicate 
if  the  top  is 
burning , 


FUNCTION  DISPLAY 


XEQ[CLRG] 


E+ 

THIS  PROGRAM 
DETERMINES  THE 
GEOMETRY  OF 

CP  GRAIN  WITH 

AN  OFF  CENTER 

OR  COCKED 

MANDREL 

R/S 

COCKED? 

Y=l,  N=0 

R/S 

R/S 

COCKED  AT  TOP  ONLY 
Y=l,  N=0 

R/S 

LGRAIN  =  ? 

R/S 

R/S 

BOTTOM  BURNING 

Y=l,  N=0 

R/S 

TOP  BURNING 

Y=l,  N=0 

R/S 

TOP  BURNING 

Y=>1,  N=0 

R/S 

LGRAIN  =  ? 

Table  A-1.  Program  Instructions  -  Continued 


STEP 

INSTRUCTIONS 

INPUT 

FUNCTION 

DISPLAY 

6. 

Key  in  grain 
length. 

L(0) 

R/S 

RGRAIN  =  ? 

7. 

Key  in  grain 
radius. 

R(0) 

R/S 

7. a 

If  grain  is 
cocked  go  to  Step  8. 

OFF  SET  =  ? 

7.b 

If  grain  is 
not  cocked. 

NO.  END  BURN  »  ? 

7.b.l 

If  grain  is 
not  cocked  enter 
number  of  ends 
burning* 

Neb 

R/S 

OFF  SET  =  ? 

8. 

Key  in  mandrel 
off  set. 

AX  or  AXt 

R/S 

TAU  START  =  ? 

9. 

Key  in  starting 
web  distance 
burned. 

■'^start 

R/S 

TAU  STOP  =  ? 

10. 

Key  in  stopping 
web  distance 
burned. 

'’^stop 

R/S 

10. a 

If  start  and 
stop  are  equal 
go  to  step  11. 

10. b 

If  start  and 
stop  are  not 
equal. 

DELTA  TAU  %  =  ? 

lO.b.l 

If  start  and 

At 

R/S 

stop  are  not 
equal  enter  the 
web  distance  increment 
then  go  to  step  11. 

11.  Write  program 

run  Information. 

11. a  If  mandrel  is  not  GEOMETRY  FOR 

cocked.  CP  GRAIN  WITH 

AN  OFFSET  OF 
X.XXXXX  IN 
AND  X  ENDS  BURNING 
SHORT  WEB  =•  X.XXXXXX 
MAX  WEB  -  X.XXXXXX 

A-1 


F 

i 


U 


STEP 


Table  A-1.  Program  Instructions  -  Continued  j 


INSTRUCTIONS _ INPUT _ FUNCTION  DISPLAY 


11. b 

If  mandrel  is 

GEOMETRY  FOR  COCKED 

cocked  at  the 

CP  GRAIN  WITH 

top  and  bottom. 

AN  OFF  SET  OF 

X.XXXXX  IN 

AND  X  ENDS  BURNING 
SHORT  WEB  =  X.XXXXXX 
MAX  WEB  =  X.XXXXXX 

11. c 

If  mandrel  is 

GEOMETRY  FOR  COCKED 

cocked  at  the 

AT  TOP  ONLY 

top  only. 

CP  GRAIN  WITH 

AN  OFFSET  OF 

X.XXXXXX  IN 

AND  X  ENDS  BURNING 
SHORT  WEB  =  X.XXXXXX 
MAX  WEB  =  X.XXXXXX 

12. 

Display  motor 

TAU  =  X.XXXXXX  IN 

geometries 

%  Web  =  XX.XXXX% 

for  web  distance 

Ab  =  XXX.XXXX  SQ  In 

burned  values  from 

"^start  ■’■stop 

in  Increments  of 

At.  Also  display 

VOL  =  XXX.XXXX  CU  IN 

the  geometry  for 

TAU  =  X.XXXXXX  IN 

the  point  of  short 

%  Web  =  XX.XXXX% 

web  burn  out.  In 

Ab  =  XXX.XXXX  SQ  IN 

addition  program 
will  stop  at  t„5q 

VOL  =  XXX.XXXX  CU  IN 

if  tstop  exceeds 

SHORT  WEB 

%bo. 

BURN  OUT 

13. 

(Optional) 

Evaluate  a 
single  motor 
geometry. 

1/X  TAU  =  ? 

14. 

(Optional)  Key 

R/S  TAU  =  X.XXXXXX  IN 

in  web  distance 

%  WEB  =  XX.XXXX% 

burned  to  be 

Ab  =  XXX.XXXX  SQ  IN 

evaluated , 

VOL  =  XXX.XXXX  CU  IN 

15. 

(Optional) 

Evaluate  the  same 
motor  geometry  with 
a  new  offset  value. 

Return  to  step  8. 

OFF  SET  =  ? 
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To  evaluate  a  new 
problem  go  to  step  3 


When  operating  the  program,  note  that  as  long  as  the  program  registers  are 
not  cleared  all  Input  values  are  maintained  until  they  are  specifically 
replaced.  If  any  portion  of  the  input  sequence  is  initiated,  the  previous  value 
for  any  input  variable  will  be  retained  if  R/s  is  entered  after  the  respective 
prompt.  Thus,  for  an  input  value  to  be  changed  at  a  prompt,  a  numeric  entry 
must  be  made. 

Another  item  that  should  be  noted  when  operating  the  program  is  the  value  of 
At.  If  the  mandrel  is  not  cocked  the  value  of  Ax  is  the  input  as  a  percentage 
of  Xpbo*  If  the  mandrel  is  cocked.  Ax  is  input  as  a  percentage  of  Xmbo* 
addition,  if  the  mandrel  is  not  cocked,  the  short  web  value  that  is  output  is 
Xg,,  and  the  maximum  web  value  is  Xpbo*  If  the  mandrel  is  cocked  the  short  web 
value  is  xj  and  the  maximum  web  value  is  x^bo* 

B.  Sample  Problems 

With  the  operation  of  the  program  completely  detailed,  the  next  step  is  to 
demonstrate  the  use  of  the  program  on  some  sample  problems.  For  a  sample  motor 
geometry  the  2X4  ballistic  test  motor  was  chosen.  The  basic  dimensions  of 
this  motor  are  as  follows: 

Rc  =  1.00  in. 

L(0)  =  3.75  in. 

R(0)  =  .75  in. 

In  the  first  sample  problem,  the  program  exercised  was  for  a  grain  configuration 
which  was  cast  with  a  mandrel  cocked  at  the  top  only.  For  this  same  geometry 
the  "ONE"  and  "START"  options  were  also  demonstrated.  The  "START"  program 
option  allows  the  user  to  evaluate  the  same  basic  configuration  with  a  different 
degree  of  mandrel  misalignment.  The  "ONE”  option  allows  the  user  to  evaluate 
the  present  configuration  at  a  single  web  distance  burned.  The  program  was  also 
exercised  for  two  other  grain  configurations,  a  grain  cast  with  a  mandrel  cocked 
at  both  the  top  and  bottom  and  a  grain  cast  with  a  displaced  mandrel.  The 
complete  details  of  these  sample  problems  are  as  follows: 


Mandrel  Cocked  at  Top  Only 


XEO  'OFCHTR- 

THIS  PROGRflH 
DETERHIHES  THE 
GEOHETRY  OF 
CP  GRRIM  UITH 
m  OFF  CENTER 
OR  COCKEI 
WMBREL 

RCRSE='> 

l.MMMMe  RUN 

COCKEI’ 

Y=l,  H=0 

QI  ly 

1  •  WWWWW  Kur 

COCKEI  AT  TOP  ONLY 
Y=li  H=8 

QHhi 

•  RUn 

BOTTOM  BURNING’ 

Y=l.  N=B 

l.MBMBBM  RUN 
TOP  BURNING’ 

Y=l,  H=0 

1.000000000  RUN 

LGR(IIN=’ 

3.7S0000000  RUN 

RGRRIN-? 

TCiMAAAAA  Dliy 

•  f  jWIPVvWv  Run 

OFF  SET=’ 

.040000000  RUN 

TRU  8TfiRT=’ 

.200000000  RUN 

TflU  ST0P=? 

.220000000  RUN 

BELTR  TRU  •/.=? 

2.000000000  RUN 

GEOMETRY  FOR 
COCKED 
RT  TOP  ONLY 
CP  GRAIN  HITH 
RN  OFFSET  OF  0.04000  IN 
AND  2.  ENDS  BURNING 
SH(»T  HEB=0.212264  IN 
MAX  HEB=0.286939  IN 


TRU=0. 200000  IN 
y.  HEB=69.7011  V. 
Rb=20.6088  SQ  IN 
VOL=10.7548  CU  IN 


TflO=0. 205739  IN 
y.  HEB=7i.70ii  y. 
Rb=20.5920  SQ  IN 
VOL=10.8731  CU  IN 

TfiU=0.211478  IN 
y.  HEB=73.70U  7, 
Rb=20.5739  SQ  IN 
VOL=10.9912  CU  IN 

TRO=0.2I2264  IN 
7.  HEB=?3.9753  y. 
Rb=28.57n  SQ  IN 
VOL=11.0074  CU  IN 

snort  neb 

BURN  OUT 


TRU=0.217216  in 
>.  HEB=75.70U  7. 
ftb=20.2429  SO  IN 
VOl=11.1088  CU  IH 

TRU=0. 220000  IN 
7.  HEB=76.8713  7. 
«b=19.922l  SQ  IN 
V0L=l 1.1645  CU  IN 


"START"  for  Mandrel 
Cocked  at  Top  Only 

XEQ  ’STRRT’ 

OFF  SET=’ 

.035000000  RUN 

TflU  STflRT=? 

.230000000  RUN 

TflU  ST0P=’ 

.240000000  RUN 
lELTft  TflU  7.=? 

5.000000000  RUN 

GEOMETRY  FOR 
COCKED 
AT  TOP  ONLY 
CP  GRAIN  WITH 
AN  OFFSET  OF  0.03500  IN 
AND  2,  ENDS  BURNING 
SHORT  HEB=0.217026  IN 
MAX  HEB=0. 282365  IN 
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TnU=0. 230866  IN 
%  UEB=81.45S0  •/. 
«b=18.7354  SQ  IN 
V0L=1 1.3627  CU  IN 

TflU=0. 240060  IN 
HEB=84.9%5  ’/. 
ftb=15.6999  SQ  IN 
V0L=11.5358  Cl)  IN 


"ONE"  for  Mandrel 
Cocked  at  Top  Only 

XEQ  -ONE" 

TflU=? 

.260008600  RUN 


TflU=0. 266000  IN 
■/.  HEB=92.0795  ’4 
nb=4.7974  SQ  IN 
V0L=1 1.7406  CU  IN 


Mandrel  Cocked  at 
Top  and  Bottom 


XEQ  -OFCNTR- 

THIS  PROGRQH 
BETERHINES  THE 
GEOHETRY  OF 
CP  GROIN  HITH 
ON  OFF  CENTER 
OR  COCKED 
nONDREL 

RCflSE=? 

1.000000000  RUN 

COCKED? 

Y=l,  N=0 

1.000000000  RUN 
COCKED  OT  TOP  ONLY 
Y=l,  N=0 

0.000000000  RUN 
TOP  BURNING? 

Y=l,  N=0 

1.000000000  RUN 
LGR0IH=? 

3.750800000  RUN 
RGRflIN=? 

.750000000  RUN 


OFF  SET=? 

.040600000  RUN 
TflU  STORT=? 

.200000000  RUN 
TOU  ST0P=? 

.220000000  RUN 
DELTA  TOU  ?.=? 

2.000000000  RUN 

GEOMETRY  FOR 
COCKED 

CP  GROIN  WITH 
ON  OFFSET  OF  0.04006  IN 
AND  2.  ENDS  BURNING 
SHORT  UEB=0.214S78  IN 
MAX  HEB=0. 283943  IN 


TOU=0. 200000  IN 
y.  UEB=70.4368  >. 
Ob=20.6088  SQ  IN 
VOL=10.7548  CU  IN 

TflU=0. 205679  IN 
::  HEB=72.4368 
Ob=20.5922  SQ  IN 
VOL=10.8718  CU  IN 

TOU=0.211358  IN 
'/.  NEB=74.4368  5; 
Ob^20.5743  SQ  IN 
VOL=10.9887  CU  II' 

TOU=0. 214578  IN 
■/.  NEB=75.5708  >. 
Ob=20.5636  SQ  IN 
VOL=11.0550  CU  IN 

SHORT  NEB 
BURN  OUT 


TOO=0.217037  IN 
"/.  UEB=76.4368 
Ob=20.4475  SQ  IN 
VOL=11.1054  CU  IN 

THU=0. 220000  IN 
y.  HEB=77.4805  !! 
Ob=20.1830  SQ  IN 
V0L=1 1.1656  CU  IN 
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Displaced  Mandrel 


SHORT  HEB 
BURN  OUT 


THIS  PROGRAM 

TAU=0.211600  IN 

DETERMINES  THE 

X  HEB=72.9655  X 

GEOMETRY  OF 

Ab=18.7221  SQ  IN 

CP  GRAIN  HITH 

VOL=10.9917  CU  IN 

AN  OFF  CENTER 

OR  COCKED 

TAU=e.217400  IN 

MANDREL 

HEB=74.9655  X 

Ab=lB.5252  SQ  IN 

RCASE-? 

V0L=1 1.0932  CU  IN 

1.000008000 

RUN 

COCKED? 

TRU=0. 220000  IN 

Y=l.  M=0 

X  HEB=75.8621  X 

0.080000000 

RUN 

Ab=15.8357  SQ  IN 

LGRAIN='> 

V0L=11.1353  CU  IN 

3.750000000 

RUN 

RGRAIN=? 

.750000000 

RUN 

HO.  END  BURH=? 

2.000000000 

RUN 

OFF  SET=? 

.040000000 

RUN 

TAU  START=? 

.200000000 

RUN 

TAU  STOP=? 

.220000800 

RUN 

DELTA  TflU  •/.=■? 

2.000080000  RUN 

GEOMETRY  FOR 
CP  GRAIN  WITH 
AN  OFFSET  OF  0.04088  IN 
AND  2.  ENDS  BURNING 
SHORT  UEB=O.210OOO  IN 
MAX  UEB=0. 290000  IN 


TAU=8. 200000  IN 
'/.  HEB=68.9655  X 
Ab=20.6088  SQ  IN 
VOL=10.7548  CU  IN 

TAU=0. 205800  IN 
X  HEB=70.9655  X 
Rb=20.5918  SQ  IN 
VOL=10.8743  CU  IN 

TAU=0. 210000  IN 
X  HEB=72.4138  X 
Ab=20.5787  SC  IN 
VOL=10.9608  CU  IN 
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C.  Installation  Information 


With  the  operational  aspects  of  the  program  presented,  the  next  step  Is  to 
provide  the  information  required  to  install  the  program  on  an  HP-41C  calculator 
system.  Presented  below  is  a  complete  listing  of  the  program.  From  this 
listing  the  program  can  be  directly  keyed  into  the  calculator.  To  facilitate  an 
understanding  of  the  program  listing,  the  storage  reslster  assignments  are  pre¬ 
sented  in  Table  A-2,  and  to  aid  in  the  installation  and  operation  of  the  program 
information  about  the  required  calculator  status  is  presented  In  Table  A-3. 


ei^LBL  "OFC 
NTR" 

02  FIX  9 
03  CF  22 
04  CF  01 
05  CF  02 
06  CF  03 
07  0DV 


08 

"THIS  PR 

OGROM" 

09 

HVIEW 

10 

"DETERM  I 

NES 

THE" 

1  1 

PVIEW 

12 

"GEOMETR 

Y  OF" 

13 

PVIEW 

14 

"CP  GRP  I 

N  WITH" 

15 

PVIEW 

16 

"PH  OFF 

CENTER" 

17 

PVIEW 

18 

"OR  COCK 

ED" 

19 

PVIEW 

20 

"MPNDREL 

21 

PVIEW 

22 

PDV 

23 

"RCPSE=? 

24 

PROMPT 

25 

FS?  22 

26 

STO  09 

27 

CF  22 

28 

"COCKED? 

29 

PVIEW 

30 

"Y=l>  N= 

Input 


0 


31 

PROMPT 

32 

FS?  22 

33 

STO  40 

34 

CF  22 

35 

1 

36 

RCL  40 

37 

- 

38 

X<=0? 

39 

SF  02 

40 

FC?  02 

41 

GTO  76 

42 

"COCKED 

PT  TOP  0" 

43 

"l-HLY" 

44 

PVIEW 

45 

"Y=l.  N 

0" 

46 

PROMPT 

47 

FS?  22 

48 

STO  42 

49 

CF  22 

50 

1 

51 

RCL  42 

52 

- 

53 

X<=0? 

54 

SF  03 

55 

0 

56 

FC?  03 

57 

STO  37 

58 

FC?  03 

59 

GTO  65 

60 

"BOTTOM 

BURNING?" 

61 

PVIEW 

62 

"Y=l.  N 

0" 

63 

PROMPT 

64 

FS?  22 

A-9 


65 

STO  37 

115 

X<=0? 

66 

CF  22 

116 

GTO  1  1 

67^LBL  65 

117 

0 

68 

■TOP  BUR 

118 

STO  17 

Evaluate  the 

NING?  ■ 

119^LBL  11 

same  motor 

69 

nviEUi 

120^LBL  "STO 

conf iguration 

70 

••Y=l.  N= 

RT- 

with  a  new 

0" 

Input 

121 

CF  01 

mandrel  offset 

71 

PROMPT 

122 

•OFF  SET 

72 

FS?  22 

=  *?  ” 

73 

STO  36 

123 

PROMPT 

74 

CF  22 

124 

FS?  22 

75^LBL  76 

125 

STO  38 

76 

■■LGRfiIN  = 

126 

CF  22 

? 

127 

RCL  38 

77 

PROMPT 

128 

STO  08 

78 

FS?  22 

129 

"TOU  STO 

79 

STO  02 

RT^'T 

*  •« 

80 

FC?  22 

130 

PROMPT 

81 

GTO  43 

131 

FS?  22 

82 

FC?  02 

132 

STO  21 

83 

GTO  43 

133 

CF  22 

84 

2 

134 

•TOU  STO 

85 

FC?  03 

P=?- 

Input 

86 

ST/  02 

135 

PROMPT 

87^LBL  43 

136 

FS?  22 

88 

CF  22 

137 

STO  22 

89 

■■RGRRIN  = 

138 

CF  22 

?  ” 

139 

0 

90 

PROMPT 

140 

STO  00 

91 

FS?  22 

141 

XEQ  "GEO 

92 

STO  01 

•a 

93 

CF  22 

142 

FS?  02 

94 

FS?  02 

143 

XEQ  "GEO 

95 

GTO  1  1 

2" 

96 

••NO.  END 

144 

RCL  22 

BURN=?" 

145 

RCL  21 

97 

PROMPT 

146 

98 

FS?  22 

147 

X<=0? 

99 

STO  17 

148 

GTO  20 

100 

CF  22 

149 

"DELTO  T 

101 

RCL  17 

ou  : 

102 

INT 

150 

PROMPT 

103 

STO  17 

151 

FS?  22 

104 

2 

152 

STO  20 

105 

RCL  17 

153 

FC?  22 

106 

— 

154 

GTO  20 

107 

CHS 

155 

RCL  20 

108 

X<=0? 

156 

100 

109 

GTO  10 

157 

110 

2 

158 

RCL  16 

111 

STO  17 

159 

* 

112^LBL  10 

160 

FC?  02 

Input 

113 

RCL  1  7 

161 

GTO  89 

114 

CHS 

162 

RCL  16 

A-10 


163  / 

164  RCL  27 

165  ♦ 

166^LBL  89 
167  STO  20 
168^LBL  20 

169  CF  22 

170  ODV 

171  "GEOMETR 
Y  FOR" 

172  OVIEW 

173  FC?  02 

174  GTO  79 

175  "COCKED" 

176  fiVIEW 

177  FC?  03 

178  GTO  79 

179  “OT  TOP 
ONLY" 

180  nviEU 
ISl^LBL  79 

182  "CP  GRfil 
N  WITH" 

183  fiVIEW 

184  FIX  5 

185  RCL  08 

186  FS?  02 

187  RCL  38 

188  "ON  OFFS 
ET  OF  " 

189  ORCL  X 

190  “F  in- 

191  OVIEW 

192  FIX  0 

193  RCL  17 

194  FC?  02 

195  GTO  80 

196  2 

197  RCL  36 

198  RCL  37 

199  + 

200  FC?  03 

201  * 

202^LBL  80 

203  "OND  " 

204  ORCL  X 

205  "F  ENDS 
BURNING” 

206  OVIEW 

207  FIX  6 

208  RCL  10 

209  FS?  02 

210  RCL  24 

211  "SHORT  W 


Output 
information 
on  the  case 
being  evaluated 


EB=" 

212  ORCL  X 

213  "  F  I  N 

214  OVIEW 

215  RCL  16 

216  FS?  02 

217  RCL  27 

218  “MOX  WEB 

219  ORCL  X 

220  "F  IN- 

221  OVIEW 
222^LBL  95 

223  FIX  9 

224  ODV 

225  RCL  21 

226  STO  00 

227  FC?  02 

228  RCL  10 

229  FS?  02 

230  RCL  24 

231  - 

232  CHS 

233  X<=0? 

234  SF  01 
235^LBL  30 

236  FS?  02 

237  XEQ  "GEO 

•• 

238  FS?  02 

239  GTO  81 

240  XEQ  "GEO 

•• 

241^LBL  81 

242  XEQ  "OUT 
PUT" 

243  RCL  22 

244  RCL  21 

245  - 

246  X<=0? 

247  STOP 

248  RCL  20 

249  ST+  00 

250  RCL  22 

251  RCL  00 

252  - 

253  X<=0? 

254  GTO  35 

255  RCL  16 

256  FS?  02 

257  RCL  27 

258  RCL  00 

259  - 

260  X<=0? 


261 

GTO 

31 

262 

FS? 

oi 

263 

GTO 

30 

264 

RCL 

lO 

265 

FS? 

02 

266 

RCL 

24 

267 

RCL 

00 

268 

— 

269 

X<=0? 

270 

GTO 

32 

271 

GTO 

30 

272^LBL 

31 

273 

RCL 

16 

274 

FS? 

02 

275 

RCL 

27 

276 

STO 

0O 

277 

FS? 

02 

278 

2“ 

XEQ 

-GEO 

279 

FS? 

02 

280 

GTO 

82 

281 

XEQ 

-GEO 

282^LBL 

82 

283 

XEQ 

-OUT 

PUT- 

284 

GTO 

75 

285^LBL 

32 

286 

RCL 

00 

287 

STO 

23 

288 

SF  01 

289 

RCL 

10 

290 

FS? 

02 

291 

RCL 

24 

292 

RCL 

00 

293 

— 

294 

X  =  0? 

295 

GTO 

47 

296 

RCL 

10 

297 

FS? 

02 

298 

RCL 

24 

299 

STO 

00 

300 

FS? 

02 

301 

2- 

XEQ 

-GEO 

302 

FS? 

02 

303 

GTO 

83 

304 

XEQ 

-GEO 

305^LBL 

83 

306 

XEQ 

"OUT 

PUT 

~ 

307^LBL 

47 

308 

RCL 

23 

309 

STO 

00 

TIC 

QDV 

31  1 

-SHORT  U 

EB- 

312 

QVIEU 

313 

"BURN  OU 

T“ 

314 

nviEM 

315 

QDV 

316 

GTO 

30 

317^LBL 

35 

318 

RCL 

22 

319 

STO 

00 

320 

FS? 

02 

321 

XEQ 

-GEO 

2- 

322 

FS? 

02 

323 

GTO 

84 

324 

XEQ 

-GEO 

325^LBL 

84 

326 

XEQ 

-OUT 

PUT- 

327 

STOP  - 

328^LBL 

-OUT 

PUT- 

329 

1 

330 

FC? 

03 

331 

2 

332 

FC? 

02 

333 

1 

334 

ST* 

13 

335 

ST* 

14 

336 

QDV 

337 

FIX 

g  Output  the  motor 

338 

RCL 

00  geometry  for  the 

339 

-  given  web  die- 

340 

QRCL  X  tance  burned 

341 

-H 

IN- 

342 

QVIEU 

343 

RCL 

OO 

344 

FC? 

02 

345 

RCL 

16 

346 

FS? 

02 

347 

RCL 

27 

348 

✓ 

349 

lOO 

350 

351 

FIX 

4 

352 

WEB=" 

353 

QRCL  X 

354 

-h 

5:- 

355 

RVIEM 

356 

RCL 

13 

357 

-fib 

358 

fiRCL  X 

359 

"h 

SQ  IN 

A-12 


360  0VIEU 

361  RCL  14 

362  ••VOL=" 

363  RRCL  X 

364  "I-  CU  IN 

365  RVIEU 

366  FIX  9 

367  RTN 
368^LBL  "GEO 

369  RCL  01 

370  RCL  00 

371  + 

372  STO  03 

373  RCL  02 

374  RCL  00 

375  RCL  17 

376  * 

377  - 

378  FC?  02 

379  STO  04 

380  RCL  09 

381  RCL  08 

382  + 

383  RCL  01 
304 

385  STO  16 

386  RCL  16 

387  FS?  02 

388  RCL  27 

389  RCL  00 

390  - 

391  CHS 

392  X<=0? 

393  GTO  01 
394^LBL  75 

395  ODV 

396  "MOTOR  B 
URNED" 

397  RVIEM 

398  "OUT" 

399  RVIEM 

400  PDV 

401  STOP 
402^LBL  01 

403  RCL  09 

404  RCL  01 

405  - 

406  RCL  08 

407  - 

408  STO  10 

409  RCL  00 

410  - 


Calculate  the 
geometry  for 
a  motor  cast 
with  an  offset 
mandrel 


41  1 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 


X<0? 
GTO  00 
RCL  09 
XT2 

RCL  03 
Xt2 


RCL  17 
* 

RCL  03 
2 


RCL  04 


STO  13 
RCL  09 
XT2 
PI 
♦ 

RCL  02 
* 

RCL  09 
XT2 

RCL  03 
Xt2 


RCL  04 


STO  14 
360 

STO  06 
STO  07 
RCL  03 
2 


STO  05 
RCL  09 
XT2 

RCL  03 
Xt2 


STO  19 
RTN 


466^LBL 

OO 

52  1 

.X 

467 

RCL 

03 

522 

RTRN 

468 

X1‘2 

523 

2 

469 

RCL 

09 

524 

3|t 

470 

X1‘2 

525 

STO  07 

471 

— 

526 

GTO  04 

472 

RCL 

08 

527^LBL  03 

473 

xrz 

528 

RCL  18 

474 

— 

529 

RCL  15 

475 

2 

530 

y 

476 

y' 

531 

CHS 

477 

RCL 

08 

532 

RTHH 

478 

y 

533 

2 

479 

STO 

15 

534 

♦ 

480 

XT2 

535 

CHS 

481 

CHS 

536 

360 

482 

RCL 

09 

537 

+ 

483 

X1‘2 

538 

STO  07 

484 

+ 

539^LBL  04 

485 

RBS 

540 

RCL  15 

486 

SORT 

541 

RCL  08 

487 

STO 

18 

542 

+ 

488 

RCL 

08 

543 

X^0? 

489 

RCL 

03 

544 

GTO  05 

490 

+ 

545 

180 

491 

RCL 

09 

546 

STO  06 

492 

+ 

547 

GTO  07 

493 

2 

548^LBL  05 

494 

y 

549 

RCL  15 

495 

STO 

1  1 

550 

RCL  08 

496 

RCL 

03 

551 

+ 

49^7 

— 

552 

X<0? 

498 

RCL 

1 1 

553 

GTO  06 

499 

RCL 

08 

554 

RCL  08 

500 

— 

555 

RCL  15 

501 

4> 

556 

+ 

502 

RCL 

1 1 

557 

1-^X 

503 

RCL 

09 

558 

RCL  18 

504 

— 

559 

4c 

505 

* 

560 

RTRN 

506 

RCL 

1  1 

561 

2 

507 

562 

♦ 

508 

SORT 

563 

STO  06 

509 

STO 

12 

564 

GTO  07 

510 

RCL 

15 

565^LBL  06 

51 1 

X*0? 

566 

RCL  08 

512 

GTO 

02 

567 

RCL  15 

513 

180 

568 

+ 

514 

STO 

07 

569 

CHS 

515 

GTO 

04 

570 

l^X 

516^LBL 

02 

571 

RCL  18 

517 

X<0? 

572 

4c 

518 

GTO 

03 

573 

RTRN 

519 

RCL 

18 

574 

2 

520 

RCL 

15 

575 

4c 

576 

CHS 

63 1 

577 

360 

632 

CHS 

578 

+ 

633 

PI 

579 

STO 

06 

634 

RCL 

09 

58e^LBL 

07 

635 

XT2 

581 

PI 

636 

a|c 

582 

180 

637 

RCL 

02 

583 

y 

638 

a|t 

584 

RCL 

06 

639 

+ 

585 

640 

STO 

14 

586 

RCL 

03 

641 

RCL 

09 

587 

* 

642 

XT2 

588 

STO 

05 

643 

RCL 

07 

589 

RCL 

09 

644 

4c 

598 

XT2 

645 

RCL 

03 

591 

RCL 

07 

646 

XT2 

592 

«c 

647 

RCL 

06 

593 

RCL 

03 

648 

4c 

594 

XT2 

649 

— 

595 

RCL 

06 

650 

PI 

596 

* 

651 

4c 

597 

— 

652 

360 

598 

PI 

653 

y 

599 

654 

RCL 

12 

608 

360 

655 

2 

601 

y 

656 

4c 

602 

RCL 

12 

657 

+ 

603 

2 

658 

STO 

19 

604 

659 

RTN 

605 

+ 

660^LBL 

-ONE 

606 

RCL 

17 

M 

Calculate  the 

607 

4t 

661 

CF  01 

motor  geometry 

608 

RCL 

04 

662 

"TOU=?“ 

for  a  single 

609 

RCL 

05 

663 

PROMPT 

web  distance 

610 

>K 

664 

FS? 

22 

burned 

61 1 

+ 

665 

STO 

21 

612 

STO 

13 

666 

CF  22 

613 

RCL 

09 

667 

RCL 

21 

614 

XT2 

668 

STO 

22 

615 

RCL 

07 

669 

1O0 

616 

670 

STO 

20 

617 

RCL 

03 

671 

GTO 

95 

618 

XT2 

672^LBL 

-GEO 

619 

RCL 

06 

2" 

620 

a* 

673 

RCL 

09 

621 

— 

674 

RCL 

01 

622 

PI 

675 

— 

623 

676 

RCL 

38 

624 

360 

677 

— 

625 

y 

678 

1 

626 

RCL 

12 

679 

RCL 

:',s 

627 

2 

680 

RCL 

36 

628 

6S1 

629 

+ 

y  2.' 

F  C  L 

Cl 

630 

RCL 

04 

A 

79-4 

PCL 

09 

P" 

7^S 

RC'_ 

81 

848 

RTK 

796 

— 

8494LBL 

73 

797 

RCL 

00 

850 

RCL 

00 

798 

— 

851 

RCL 

27 

799 

RCL 

38 

852 

— 

8^0 

853 

X>0? 

881 

RCL 

82 

854 

GTO 

75 

882 

855 

RCL 

00 

883 

STO 

29 

856 

RCL 

09 

804 

STO 

41 

857 

— 

885 

RCL 

82 

858 

RCL 

01 

886 

RCL 

80 

859 

+ 

887 

RCL 

36 

868 

RCL 

02 

888 

861 

4 

889 

— 

862 

RCL 

38 

818 

STO 

38 

863 

✓ 

81 1 

RCL 

28 

864 

STO 

28 

812 

— 

865 

STO 

41 

813 

STO 

84 

866 

RCL 

02 

814 

XEQ 

"RRE 

867 

RCL 

00 

R“ 

868 

RCL 

36 

815 

RCL 

29 

869 

4 

816 

RCL 

28 

870 

— 

817 

— 

871 

STO 

30 

818 

STO 

23 

872 

RCL 

28 

819 

RCL 

31 

873 

— 

820 

4 

874 

STO 

04 

821 

ST  + 

13 

875 

XEQ 

“ORE 

822 

RCL 

23 

R- 

823 

RCL 

34 

876 

RTN 

824 

lie 

8774LBL 

“RRE 

825 

ST- 

14 

R“ 

826 

RTN 

878 

RCL 

28 

827^LBL 

72 

879 

RCL 

02 

828 

RCL 

08 

880 

✓ 

829 

RCL 

26 

881 

RCL 

38 

838 

— 

882 

4 

831 

X>0? 

883 

STO 

08 

832 

GTO 

73 

884 

XEQ 

-GEO 

833 

RCL 

37 

as 

834 

RCL 

80 

885 

RCL 

05 

835 

4 

886 

STO 

31 

836 

STO 

28 

887 

RCL 

19 

837 

STO 

41 

888 

STO 

34 

838 

RCL 

02 

889 

RCL 

30 

839 

RCL 

00 

890 

RCL 

02 

848 

RCL 

36 

891 

z' 

841 

4 

892 

RCL 

38 

842 

— 

893 

4 

843 

STO 

30 

894 

ST  J 

08 

b44 

RCL 

28 

895 

XEQ 

-GEO 

345 

— 

as 

STO 

04 

89t. 

PCL 

05 

947 

xe* 

"PPE 

89  7 

STO 

Integrate  to 
obtain  surface 
area  and  volume 
using  a  trape¬ 
zoidal  rule. 
Approximation 
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898 

RCL 

19 

935 

RCL 

19 

899 

STO 

35 

936 

ST  + 

39 

900 

RCL 

31 

937 

RCL 

17 

901 

RCL 

33 

938 

ST+ 

41 

902 

+ 

939 

CF  21 

903 

2 

940 

VIEU 

1  41 

904 

y 

941 

SF  21 

905 

STO 

32 

942 

RCL 

41 

906 

RCL 

34 

943 

1 . 00001 

907 
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TABLE  A-2.  Register  Assignments 
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deg 
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In 

12 
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13 
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Xl 
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18 

Yl 
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At 

In 

21 
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24 

•>^1 
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T2 
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26 

•^3 
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27 

■'inbo 

in 

28 
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in 

29 

^ub 

in 

30 

^top 

in 

31 

in 
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P(Zub)»  P 

in,  in 

33 

P( Zfop) 

in 

34 

35 

36 
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^op 

in 

in 

NA 

37 

^bot 

NA 

38 

AXt 

In 

39 

Z  Acr 

in 

40 

used 

NA, 

41 

Z 

in  2 

42 
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NA 
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TABLE  A-3.  Calculator  Status 


Calculator 

mode 

USER 

Size 

43 

Program 

registers 

276 

Total 

registers 

319 

n 

4J 

G 

0) 

E+ 

OFCNTR 

§ 

1/X 

ONE 

00 

< 

X 

START 

Flag 

Set  Flag 

Cleared  Flag 

No. 

Indicates 

Indicates 

00 

D 

u 

01 

Web  distance  burned 

Web  distance  burned  has 

cd 

has  exceeded  the 

not  exceeded  the  short 

t/i 

short  web 

web 

00 

<d 

u* 

02 

Mandrel  Is  Cocked 

Mandrel  Is  not  cocked 

03 

Mandrel  Is  cocked 

Mandrel  Is  cocked 

at  top  only 

top  and  bottom 
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APPENDIX  B 


MISALIGNED  2X4  MOTOR 

The  2X4  ballistic  test  motor  Is  the  basic  burning  rate  characterization 
motor  employed  by  the  Propulsion  Directorate.  This  motor  has  a  cylindrical  port 
and  In  normal  applications  has  both  end  surfaces  uninhibited. 

Initial  port  radius: 

R(0)  =■  .75  In 

Initial  grain  length: 

L(0)  =  3.75  In 

Case  Radius: 

R^  =  1 .00  In 

This  motor  was  used  as  an  example  to  demonstrate  the  application  of  the 
misaligned  motor  geometry  model.  For  this  motor  the  burning  surface  area 
histories  were  generated  for  geometries  reflecting  a  perfectly  aligned  mandrel, 
a  displaced  mandrel,  a  mandrel  cocked  at  the  top,  and  a  mandrel  cocked  at  both 
the  top  and  bottom.  For  the  three  modes  of  mandrel  misalignment,  surface  area 
histories  were  generated  for  AX^  values  of  0.00  in.,  0.01  in.,  0.02  in.,  0.03 
in.,  0.04  in.,  0.05  in.,  0.06  in.,  0.07  in.,  0.08  in.,  0.09  in.,  and  0.10  in. 

All  these  surface  area  histories  were  generated  using  the  HP-4 1C  calculator  and 
the  previously  detailed  program.  The  burning  surface  area  history  for  an 
aligned  2X4  motor  grain  is  presented  in  Figure  B-1 .  The  burning  surface  area 
histories  for  a  2  X  4  motor  cast  with  a  displaced  mandrel  is  presented  in  Figure 
B-2.  The  burning  surface  area  histories  for  a  2  X  4  motor  cast  with  a  mandrel 
cocked  at  the  top  is  presented  In  Figure  B-3.  The  burning  surface  area 
histories  for  a  2  X  4  motor  cast  with  a  mandrel  cocked  at  both  the  top  and  bot¬ 
tom  is  presented  in  Figure  B-4. 
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